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ABSTRACT 

We consider constraints on the formation of low-mass X-ray binaries containing 
neutron stars (NLMXBs) arising from the presence of soft X-ray transients among these 
systems. For a neutron star mass Mi ~ 1.4M0 at formation, we show that in short- 
period 1 — 2 d) systems driven by angular momentum loss these constraints require 
the secondary at the beginning of mass transfer to have a mass 1.3Mq ^ M2 ^ I.SMq, 
and to be significantly nuclear-evolved, provided that supernova (SN) kick velocities 
are generally small compared with the pre-SN orbital velocity. As a consequence a 
comparatively large fraction of such systems appear as soft X-ray transients even at 
short periods, as observed. Moreover the large initial secondary masses account for the 
rarity of NLMXBs at periods P ^ 3 hr. In contrast, NLMXB populations forming 
with large kick velocities would not have these properties, suggesting that the kick 
velocity is generally small compared to the pre-SN orbital velocity in a large fraction 
of systems, consistent with a recent reevaluation of pulsar proper motions. The results 
place also tight constraints on the strength of magnetic braking: if magnetic braking 
is significantly stronger than the standard form too many unevolved NLMXBs would 
form, if it is slower by only a factor ~ 4 no short-period NLMXBs would form at all in 
the absence of a kick velocity. The narrow range for AI2 found for negligible kick velocity 
implies restricted ranges near 4M0 for the helium star antecedent of the neutron star 
and near I8M0 for the original main-sequence progenitor. The pre-common envelope 
period must lie near 4 yr, and we estimate the short-period NLMXB formation rate in 
the disc of the Galaxy as ~ 2 x 10~^ yr~^- Our results show that the neutron star mass 
at short-period NLMXB formation cannot be significantly larger than 1.4M0. Systems 
with formation masses Mi ^ I.2M0 would have disrupted, so observations implying 
Ml ~ I.4M0 in some NLMXBs suggest that much of the transferred mass is lost from 
these systems. 



Subject headings: binaries: close — stars: neutron — stars: evolution — accretion, 
accretion discs — X-rays: stars 
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1. INTRODUCTION 

In a recent paper (King et al, 1996; hereafter Paper 1) we considered the disc instabihty model for 
soft X-ray transients (SXTs). Following van Paradijs (1996) we noted that X-ray irradiation of the 
disc surface tends to suppress the instability, so that SXT behaviour requires that the mass transfer 
rate is below a critical limit which is itself a rising function of orbital period P. For P 1 — 2 
d mass transfer is driven by the nuclear expansion of an evolved secondary, and proceeds at low 
enough rates that almost all of these systems appear as SXTs. For P ^ 1 — 2 d mass transfer is 
driven by angular momentum losses; in Paper 1 we showed that the resulting rates are higher than 
the critical rate, hence too high for SXT behaviour unless the secondary star is already somewhat 
nuclear-evolved before mass transfer begins. This effect is particularly marked in a neutron-star 
low-mass X-ray binary (neutron-star LMXB, or NLMXB): if we write (cf Paper 1) 1712 for the 
ratio of the secondary's mass to that of a ZAMS star filling the Roche lobe at the same orbital 
period, SXT behaviour requires m2 ^ 0.25. By contrast we found a much weaker limit 7712 <^ 0.75 
for SXT behaviour in a LMXB with a lOM© black-hole primary. These limits agree with the 
extreme mass ratios always found in SXTs. In Paper 1 we suggested that these limits also offered 
a potential explanation for the prevalence of black-hole systems among SXTs (8 out of 14 with 
known P) as compared with persistent LMXBs (1 out of 29 with known P). As NLMXBs can be 
reliably identified by the presence of X-ray bursts, we consider this problem further by investigating 
constraints on NLMXB formation. 



2. NLMXB FORMATION 

The requirement that short-period SXT secondaries should be significantly evolved is a powerful 
clue in investigating LMXB (and particularly NLMXB) formation. For example, the very fact that 
at least 4 NLMXBs out of 22 in the relevant period range are transients (0748-676, Cen X-4, 1658- 
298, Aql X-1; cf. van Paradijs 1995) must imply that NLMXB secondaries usually have masses 
-^2 ^ at the onset of mass transfer. Lower initial-mass secondaries would be essentially 

unevolvcd when mass transfer began, and a large population of them would require the 4 known 
SXTs among NLMXBs to be accompanied by a far greater number of persistent NLMXBs than 
observed. However to get any NLMXBs at all at P ^ 1 — 2 d simultaneously requires that some 
of the secondaries initially have masses M2 ^ 1.5M0 so that angular momentum loss via magnetic 
braking can shrink the binary. In addition, for M2 ^ 1.5M0 these NMLXBs with a main-sequence 
donor will be unstable to thermal-timescale mass transfer once contact is reached; they transfer 
mass at super-Eddington rates and do not appear as NLMXBs (Kalogera & Webbink 1996a). 
Given these restrictions all NLMXBs with P ^ 1 — 2 d (and indeed most at longer periods too) 
must apparently form with M2 in a very restricted range IMq <^ M2 <^ I.SMq. 

The properties and the very existence of this group depend sensitively on the efficiency of magnetic 
braking. The paper by Iben et al. (1995b) on general LMXB formation effectively rules out these 
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systems because it assumes that magnetic braking does not operate for M2 > IM0. Similarly 
Terman et al. (1996) do not consider them in detail as they somewhat artificially assume that 
M2 = I.OMq separates systems captured by angular momentum losses from systems evolving 
towards longer orbital period. 

Below we investigate in detail constraints on NLMXBs forming after a common envelope (CE) 
phase and subsequent helium star supernova (SN). We consider first the fundamental case of a 
spherically symmetrical SN (Sect. 2.1), and we shall show that in this case secondary stars in 
NLMXBs are indeed significantly nuclear-evolved at the onset of mass transfer. A generalization 
to the more realistic case when the neutron star receives a kick velocity during the SN (Sect. 2.2) 
reveals that the fraction of systems with unevolved secondaries increases with increasing magnitude 
of the kick. For a given magnetic braking strength this places limits on the magnitude of the kick 
velocity. Possible alternative evolutionary channels leading to NLMXB formation are discussed in 
Sect. 2.3. 

2.1. NLMXBs from spherically symmetrical helium— star supernovae 

To discover any limits on the initial mass M2 of the secondary we consider the constraints on 
the formation of LMXBs from helium-star + main-sequence (MS) binary remnants of common 
envelope evolution in the case of a spherically symmetrical SN. 

2.1.1. Constraints on the formation 

Webbink & Kalogera (1994; hereafter WK) and Kalogera & Webbink (1996a) list a number of 
requirements that the various progenitor stages must satisfy in order to lead to LMXB formation. 
Of these, three turn out to be crucial in fixing M2: 

(a) The post-common-envelope binary must be wide enough to allow the helium star to evolve to 
core collapse (requirement 6 of WK) . 

(b) The binary must survive the supernova event resulting from the core collapse of the helium star 
(7 of WK). 

(c) After the helium star progenitor of the LMXB primary explodes as a supernova the binary must 
be able to reach interaction within the age of the Galaxy (WK 8). 

If M2 is large enough and the post-SN separation is wide enough (c) will occur through the nuclear 
expansion of the secondary. Such systems require M2 ^ IM© to give main-sequence lifetimes ^ms 
less than the age tcal of the Galaxy and must have P 1 — 2 d at the onset of mass transfer; P 
will increase further as the secondary expands. However if M2 is less than this and/or the post-SN 
binary is narrower, (c) can come about through orbital angular momentum loss. After reaching 
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contact with an initial orbital period of order 1 - 2 d, such systems evolve either to longer (P 1 d) 
or shorter (P ^ 1 d) periods depending on the competition between nuclear evolution and angular 
momentum loss (Pylyser &: Savonije 1988a, b). 

For the remainder of this paper we concentrate on this latter group of systems, which are the only 
ones which can populate the short-period (P ^ 1 — 2 d) region. It is convenient to apply the 
constraints (a - c) in a different order. From (c) we require tMBj^GR < ^Gah where 



^MB = 2.2 X 10« , ^"^l,^ 



{nil + 'm2)^^^Tn'^ 



2 



yr (1) 



p \ 8/3' 



yr (2) 



77111712 

are the timescales for detached orbital evolution under magnetic braking and gravitational radiation, 
with 7 a dimensionless efficiency parameter for the former case. These shrink the post-SN binary 
from its initial period P to the value Pc at which it first comes into contact with its Roche lobe {Pd = 
P/d,nii = Ml /Mq, 7712 = M2/M0 etc). If the secondary is unevolved an adequate approximation 
for Pc is Pc — 0.375 m2 d, while Pc is up to an factor ~ 2.5 larger if the secondary is fairly massive 
(M2 ^ IMq) and close to the terminal main-sequence. 

We assume that the secondary has a structure such that a magnetic stellar wind brakes its rotation; 
since we are considering short-period LMXBs in which the secondaries are not too far from the 
main sequence this requires O.SMq ^ M2 ^ 1.5Mq. Furthermore we assume that the secondary is 
tidally locked to the orbit such that magnetic braking removes orbital (rather than only rotational) 
angular momentum. Applying simple scalings for the tidal synchronization timescale (e.g. Tassoul 
1995) suggests that this is the case for the detached systems in the period range under consideration 
(see below). The efficiency parameter 7 in (|^) allows us to test the sensitivity of our results to the 
strength of magnetic braking. The standard case 7 = 1 corresponds to the description by Verbunt 
&; Zwaan (1981) when the radius of gyration is set to \/0.2 and the calibration parameter to unity. 

Requirement (c) thus gives the longest post-SN period compatible with short-period LMXB for- 
mation. It is easy to show that the longest possible value of this period is given by equating tMB 
to tcai with P ^ Pc, so that the term in square brackets in (||) is unity. Kepler's law now gives a 
condition on the separation OpostSN of the tidally-circularized post-SN binary, i.e. 



flpostSN < 9.OP0 

where to = tGai/10^° yr. 



to (mi + 7712) ml 
7 mi 



1/5 



(3) 



Before exploiting (^) further we apply requirement (b). In a spherically symmetrical supernova 
explosion, it is well known that the remnant binary is disrupted if more than one-half of the binary 
mass is removed, i.e. unless 

MHe-Mi < i(MHe + M2) (4) 
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or 

1712 > mHe — 2mi (5) 

the binary will be disrupted. Here m-He is the helium star mass Mhc in Mq. Moreover, the 
eccentricity e of the orbit immediately after the SN is given by 

e = . (6 

mi + 7712 

Obviously condition (^) is equivalent to e < 1. The pre~SN separation and the post-SN separation 
of the tidally circularized orbit are related by (e.g. Verbunt 1993) 

OpreSN = flpostSN- (7) 

1 + e 

Using (^) and (^) we find that the pre-SN Roche lobe radius Rr of the He star must obey 



Rr ^ < 4.5Re 



1 + e 



7 m-i 



(8) 



where we have assumed that Rr occupies one-half of the pre-SN binary separation 



Finally, requirement (a) demands that the largest radius i2max(He) of the helium star must fit 
inside the pre-SN binary. The evolution of the helium star is very rapid compared with the orbital 
evolution so we assume that CpreSN is constant, and the condition is i?inax(He) < Rr. (If this 
condition fails it is possible that some LMXBs nevertheless form after a second common-envelope 
phase, cf Iben et al 1995b. We note, however, that these authors assume an unusually efficient 
envelope ejection compared to the standard common envelope description (^) below. Application 
of the latter would almost always lead to a merger. Hence we neglect this group of systems.) Iben 
et al. (1995b) fit i?max(He) as 

i?max = 77OOi?0mi^5-5 (9) 

for TTiHe > 2.5Mo. 

Hence combining (P) with (S) we see that (a) and (c) demand 



'V / o Q'7 \ 55/2 

(mi + m2)'mt > (1 + e)M ^ mi, (10) 

tc \mHeJ 



thereby defining a lower limit for m2 for given mi, mue (or a lower limit for mne for given mi, m2). 



2.1.2. Mass limits for NLMXBs 

Observational evidence strongly suggests that neutron stars emerge from the supernova that creates 
them with a mass mi = 1.4. Figure 1 shows the constraints (|5|, ^) for this mass, with = 1 and 
7 = 1. (For m2 < 0.3, where magnetic braking is thought not to operate, ( p!o| ) was replaced by the 
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corresponding constraint using tcR instead of tiviB-) As can be seen the two constraints cross at a 
value 7712 = 1-22, mHe = 4.02. Hence m2 > 1.22 for NLMXB formation through this channel. We 
now go one step further and ask under what conditions the binary can begin mass transfer with 
a secondary which has spent less than a fraction / (0 < / < 1) of its lifetime tMS on the main 
sequence, i.e. with 

^MB < / ^MS- (11) 

Clearly this is a more stringent requirement than (c), where the rhs is tcai > ^mSj St-nd it therefore 
requires a larger mass m2- Taking 

tMs = lo^Vg ^ yr (12) 

and following the same line of reasoning that led to the constraint ([lO|), i.e. replacing to in ( p!o[ ) by 

(//rn-l), we find 

'Y / 3 87 \ "'^/^ 

{mi+m2fm2>-^{l+ef(^] mi, (13) 



which must hold simultaneously with (|lO|). Condition (13) is best understood as a lower limit 
for the He star mass for any given mi, m2- In Fig. 1 we show the constraint for / = 1 (dotted 
curve) and / = 0.25 (dash-dotted), with standard magnetic braking (7 = 1). Instead of the rough 
approximation ( p^ we used the fitting formulae given by Kalogera & Webbink (1996a) to determine 



tMS- We see that condition ( p~3[ ) with / = 1 intersects the allowed mass range 1.2 ^ m2 ^ 1-5 set 
by the Galactic age constraint and the requirement of thermally stable mass transfer respectively. 
Systems with 1.2 ^ 7712 ^ 1.5 to the left of the dotted line in Fig. 1 (hatched, wide spacing) turn- 
on mass transfer with a giant donor at long orbital period (P ^ 1 — 2 d). As in most of them 
mass transfer is unstable (e.g. Kalogera & Webbink 1996a) they would probably not appear as 
long-period NLMXBs; we do not consider this group further. Systems with 1.3 ^ m2 ^ 1-5 to the 
right of this line (hatched, narrow spacing) will evolve to short orbital periods. The location of the 
dash-dotted line in Fig. 1 clearly demonstrates that all secondaries of systems in this group are 
somewhat evolved. In other words: in the limit of negligible kick velocities secondaries in a large 
fraction of NLMXBs with P ^ 1 — 2 d are close to the end of their MS lifetimes and are partially 
nuclear-evolved (rfi2 < 1) when mass transfer starts. The resulting low mass transfer rates (Pylyser 
& Savonije 1988b) then probably explain why a comparatively large fraction of NLMXBs of both 
types appear as transients. The simple scalings of Paper 1 suggest that SXT behaviour requires 
1712 ^ 0.25; full population synthesis calculations are needed in order to quantify the resulting SXT 
fraction. 

We note that Kalogera & Webbink (1996b) find that short-period LMXBs do not form at all 
without kick velocities. The reason for this difference from our conclusions is a modification of the 
standard form (|^) of magnetic braking such that ^mb varies like (rather than m^^), and the 
introduction of a reduction factor for secondaries more massive than IMq. Even if the reduction 
factor is ignored (which is ~ 7 for m2 = 1.5), their tMB is still longer by a factor ~ 6 for m2 = 1.5. 
This corresponds to 7 = 6 in ([l3|), or, equivalently, to / = 1/6. Inspection of Fig. 1 shows that the 
critical line corresponding to / = 1/6 would be to the right of the dash-dotted / = 1/4 line, hence 
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no parameter space would be left for short-period NLMBX formation. Physically this means that 
due to the correspondingly slower evolution the systems would not reach contact within the age of 
the Galaxy. 

This underlines the sensitivity of our result to the strength of magnetic braking for main-sequence 
stars in the critical mass range where the convective envelope disappears. To illustrate this further 
we consider the lower limit acrit we obtain for the post-SN orbital distance OpostSN — 2(1 + e)RR 
when the requirement Rr > iimax(He) is combined with (|9[) and (^): 

/ 4.3 

OpostSN > flcrit = 5.1i?0 (1 + e) — . (14) 

\m2 + 2mi/ 

Note that this limit does not depend on magnetic braking and that Ocrit is essentially determined 
by m2 (mi = 1.4, < e < 1). For NLMXBs with main-sequence donors of mass m2 to exist at all 
the detached evolution time tMB needed to shrink the orbit from acrit into contact must be shorter 
than iivis- On the other hand, to ensure that a large fraction of the secondaries are significantly 
evolved when mass transfer turns on, tyiB must be longer than a certain minimum fraction / ~ 1/4 
of tMS- This places a tight limit on the allowed strength of magnetic braking for systems with 
secondary mass m2, initial orbital distance Ocrit and neutron star mass mi = 1.4: 

0.25 ^ ^MB(acrit,M2) ^ ^ ^^^^ 

independent of the functional dependence of tuB on mass and period. 



2.1.3. Mass limits for progenitor systems 

The narrow range 

1-3^^2^1.5 (16) 
also implies a very restricted range for mne. from Fig. 1 we see that this must obey 

4.0 ^ mue ^ 4.3. (17) 

Using the analytic fit 

log mHe = 1-273 log mprog - 0.979 (18) 

from Kalogera & Webbink (1996b) to stellar models by Schaller et al. (1992) we see that the mass 
of the MS progenitor of the helium (and ultimately neutron) star must be in the range 

17.5 ^ mprog ^ 18.5. (19) 

Immediately before the CE phase the primary has evolved beyond core He exhaustion and has lost 
mass in the form of stellar winds. Interpolation of Schaller et al.'s models suggest that primaries 
from the mass range ( |T9| ) will then have masses /n-preCE in the range 

15.2 ^ mpreCE ^ 15.9. (20) 
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Wind mass loss carrying the specific a.m. of the mass-losing component increases the orbital sep- 
aration as the inverse of the total mass. Hence the MS parent binary orbit is closer by a factor 
~ 0.9 than the immediate pre-CE binary orbit. 

Simple treatments of common-envelope evolution (e.g. WK) lead to the relation 

2mpi.eCE("lpreCE " ^He + « Ar^j m2/2) 
apreCE — 7 «postCE, l^tj 

where apreCE, OpostCE are the binary separation before and after the common-envelope phase, a is 
the fraction of the orbital binding energy used to drive off the envelope, A is a weighting factor 
for the gravitational binding energy of the envelope to the core, and = r/j(mpreCE5 ''712) is the 
fractional size of the primary's Roche lobe at the start of the CE phase. We find with the mass 
limits (H, |0|) and r^j ~ 1/2, A ~ 1/2 



8 rUpreCE "T-preCE " '"T-He 250 

fflprcCE flpostCE- (22) 

a 1712 mue ot 

Since there is little orbital evolution before the helium star undergoes its supernova explosion, and 
the post-SN separation is within a factor of 2 of the pre-SN separation, we can translate this 
into limits on apreCE, or equivalently the orbital period of the progenitor binary. We get an upper 
limit from requirement (c) (cf |8|) , and a lower limit from requirement (a) , using (^) and the fact 
that m-He < 4.3. (This limit prevents NLMXB progenitors emerging from the common-envelope 
phase with the MS secondary already close to filling its Roche lobe, which is WK's condition 6; 
the requirement here gives a tighter limit. This in turn justifies our assumption above that the 
post-SN period P was much larger than the contact period Pc-) The resulting constraints are close, 
i.e. ~ 23OOii0/a and ~ ISOORq/u respectively. Together they show that the orbital distance Oprog 
and period -Pprog of the MS progenitor ~ ISM© + 1.4M0 binary must lie in the range 

1200 Rq ^ aprog/a ^ 2100 Rq (23) 

and 

3 yr ^ a^/^Pp,og ^ 7 yr. (24) 

The formation channel we have considered requires that the progenitor fills its Roche lobe after 
core-heluim exhaustion but before it would explode as a SN. For a ISM© MS star this implies 
1200 ^ Uprog/ Rq ^ 1600 or 3 ^ i^prog/yr ^ 4.6 (Kalogera & Webbink 1996b), thus reducing the 
limits (^3|), ( |2^ ) even further. If most short-period NLMXBs in the Galaxy form via helium-star 
supernovae, this restricted range of allowed inital orbits and the very narrow limits (H, m) then 
show that only an extremely small progenitor population can produce these systems, explaining 
their comparative rarity. Following Iben et al. (1995b) a very crude estimate of the Galactic short- 
period NLMXB formation rate i' can be obtained from 



1/ = 



(25) 
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where Ax denotes the allowed range of the quantity x. Eq. ( |25| ) relies on standard assumptions 
for the distribution of initial orbital parameters in newly forming ZAMS binaries, in particular on 
a flat distribution in log M2/Mpi.og- Using the limits obtained above we find ~ 2 x 10~'^yr^^. 
Replacing (25) by the initial parameter distribution of ZAMS binaries used in Kalogera & Webbink 
(1996b) gives a similar value for v. 



2.2. The role of kick velocities 

In reality supernova explosions are probably not spherically symmetrical, and give the compact 
primary remnant a kick velocity. In this case the main differences from the considerations in 2.1 
are twofold (see e.g. Kalogera 1996a). First, a suitably directed kick can unbind binaries that would 
have been stable according to (^), or keep together binaries that would have disrupted according to 
this inequality. The latter case requires a kick velocity comparable to the pre-SN orbital velocity 
of the companion, directed almost parallel to its instantaneous motion, and is therefore rather rare 
unless the kick velocity is close to an optimum value near the pre-SN orbital velocity. Second, 
relation (0) no longer holds, and OpostSN can be either smaller or larger than OpreSN (but never 
larger than 2api.eSN). These two effects not only widen the allowed region of NLMXB formation 
in the M2 — Mhb plane of Fig. 1, but in particular allow the formation of NLMXBs with almost 
unevolved secondaries. 

In the absence of a theoretical understanding of the origin, magnitude and direction of kick velocities 
the standard approach in population synthesis considerations (e.g. Brandt & Podsiadlowski 1995, 
Terman et al. 1996, Webbink & Kalogera 1996b) is to assume that the kicks are isotropic and that 
their magnitude derives from a given distribution function, characterized by a certain r.m.s. value 
0"^. As demonstrated e.g. extensively by Kalogera (1996a) the resulting probability distributions 
can be expressed in terms of the two governing dimensionless parameters, the ratio of kick velocity 
to relative orbital velocity in the pre-SN orbit, ^ = (Xk/worbiti and the ratio {3 of post-SN and pre-SN 
binary mass, (3 = (mi + m2)/(mHe + fn2)- 

Given the stochastic nature of the problem only a full population synthesis can provide a quantita- 
tive estimate of the fraction of NLMXBs with significantly nuclear-evolved secondaries at turn-on 
of mass transfer. To illustrate the main effect of kick velocities and to gain a very rough estimate 
of the maximum kicks we can tolerate and still maintain the large fraction of systems with evolved 
secondaries found for spherically symmetrical SNe, we make use of the analytic expression for the 
distribution of binaries over Uc = flpostSN /^preSN derived by Kalogera (1996a) under the assump- 
tion of a Maxwellian kick velocity distribution. In Fig. 2 we show the same limits as plotted in 
Fig. 1, but with the factor (1 -|- e) in ([lO|) and (|l3|) replaced by the approximate median value of 
Oc, 1.75, 1.25 and 0.75, for small (^ = 0.1), moderate (^ = 0.3) and large kick velocities (^ = 1.0), 
respectively. These median values depend only weakly on (3. For strong kicks the critical line (|5|) 
— which is equivalent to /3 = 0.5 — was replaced by the line [3 = 0.4 as the survival probability of 
these systems is only a factor of 2 smaller than the one for the most stably bound systems (corre- 
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spending to /3 ~ 0.75). The resulting enclosed area in the M2 — plane can be thought of as 
representative for the effective parameter space of NLMXB formation. Assuming that the area is a 
measure of the corresponding relative formation rate, Fig. 2 suggests that in the case of standard 
magnetic braking (|^) systems with unevolved secondaries still constitute only a small fraction of 
NLMXBs for ^ = 0.1, represent the majority for ^ = 0.3, and entirely dominate for ^ = 1.0. 

More efficient magnetic braking would increase the dominance of unevolved systems even further, 
whereas a favourable combination of weaker magnetic braking and a kick velocity distribution with 
small (or moderate) mean velocity could both ensure the formation of the (otherwise forbidden) 
class of short-period NLMXBs, and the predominance of nuclear-evolved main-sequence donors. 
However, a population subject to a large mean kick velocity would necessarily contain a large 
fraction of systems where the secondary is close to contact in the post-SN orbit and hence essentially 
unevolved at mass transfer turn-on — whatever the strength of magnetic braking. The reason for 
this is twofold. First, the limit (|l4|) formally allows post-SN orbital periods shorter than the contact 
period Pc for m-He ^ 4.7, and with strong kicks a large fraction of such binaries would survive the 
SN. Second, the SN-induced orbital reduction factor Uc is very small for the majority of systems. 

In view of this, the large fraction of soft X-ray transients observed among NLMXBs provides a 
strong argument not only against a magnetic braking stronger than our standard case (|^), but also 
against a mean kick velocity of order 350 — 400 kms~^ (^ ~ 1) invoked by Lyne & Lorimer (1994) 
from observed pulsar proper motions. Fig. 2 shows that kick velocities must on average be small 
compared to the pre-SN orbital velocity, probably ^ <^ 0.1. This is consistent with a more recent 
reevaluation of the initial velocity distribution of radio pulsars which confirms the existence of the 
high-velocity tail found by Lyne & Lorimer, but suggests that the distribution has its maximum 
at zero velocity, hence a smaller average value (Hansen & Phinney 1996, Hartmann 1996). 

2.3. Alternative evolutionary channels 

Short-period NLMXBs might also form from systems with initially fairly massive main-sequence 
secondaries above the limit for thermally stable mass transfer, 1.5 — 2 ^ 7712 ^ 3, provided these 
survive the initial phase of thermal timescale (hence super-Eddington) mass transfer. (Systems 
with m2 ^ 3 would probably develop a delayed dynamical instability, see Hjellming 1989). Kalogera 
&: Webbink (1996a, b) pointed out that such systems could reappear as stable NLMXBs with sub- 
Eddington mass transfer rate and donor mass 1 ^ m2 ^1.5. An analysis similar to the one in 2.1.2 
(using ^GR instead of ^mb) reveals that a significant fraction of them would emerge from the CE 
phase almost semi-detached. They would turn on mass transfer with an unevolved secondary, again 
in conflict with the observed comparatively large fraction of transients among NLMXBs. Hence 
we conclude that this channel cannot contribute significantly to the formation of short-period 
NLMBXs. 

Recently, Kalogera (1996b) has described yet another formation mechanism for LMXBs where no 
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CE phase is involved. Instead the orbital shrinkage is achieved by a suitably directed kick when 

the primary star in the wide progenitor binary explodes as a SN before it reaches its Roche lobe. 
Population synthesis models show (Kalogera 1996b) that the production of short-period LMXBs via 
this channel is altogether negligible compared to the standard He-star SN case if the kick velocities 
are large (cJk of order 300 — 400 kms~^), but might account for a formation rate comparable to the 
one derived in 2.1.2 if (Jk is close to an optimum value ~ 50 kms~^. However, in the latter case a 
relatively large fraction of LMXBs would start the X-ray phase with a very small (below O.SM©) 
donor mass, hence with essentially unevolved secondaries. Again, this suggests that the direct SN 
mechanism represents only a minor channel for the formation of short-period NLMXBs (whereas 
it might be important to produce long-period systems with P^i ^ 100 d) even if kick velocities are 
generally small. 



3. NLMXBs BELOW THE PERIOD GAP 

The difference between the period histograms of LMXBs and cataclysmic variables (CVs), in which 
the primary is a white dwarf, has long been remarked (e.g. White Sz Mason 1985, van den Heuvel 
&; van Paradijs 1988, Verbunt & van den Hcuvcl 1995, Kolb 1996). The most prominent difference 
is a total (or near-total) lack of LMXBs in the 80 < P < 120 min region below the famous 
CV period gap. A KS test reveals that there are also statistically significant differences above 
the gap. The hypothesis that the LMXB and CV samples are drawn from the same underlying 
distribution can be rejected at a confidence level > 99.99% in the period range 80 mins -2d 
(Fig. 3a), and at ~ 99.96% in the range 3 h - 2 d (Fig. 3b). This cannot be explained by selection 
effects discriminating against short orbital periods, since many of the X-ray periods were turned 
up by satellites such as EXOSAT which had a 4-day orbit. The differences conflict with the simple 
picture of CVs and LMXBs as essentially the same in terms of their secular evolution, apart from 
the substitution of a white dwarf by a neutron-star or black-hole primary. However our arguments 
above show that this simple picture is inaccurate, particularly for NLMXBs. CVs emerge from 
common-envelope evolution with a full range of secondary masses down to m2 ~ 0.1. The vast 
majority of these stars are essentially unevolved, and the post-CE separations are so small in many 
cases that the secondaries are close to their Roche lobes. The majority of CVs probably start 
mass transfer at periods below the gap (more than 67% according to King et al., 1994). None of 
these features hold for NLMXBs, as we have seen: the secondaries are confined to the narrow range 
1.3 ^ m2 ^ 1.5 initially, a large fraction of them must be significantly nuclear-evolved, the post CE 
(and SN) Roche lobes are considerably larger than the secondaries, and they start mass transfer at 
periods in the range 10 hr ^ P ^ 30 hr. 

We investigate differences between the CV and NLMXB period histogram which arise alone from 
these effects in Fig. 4. In particular, we test if the population of NLMXBs below P ~ 2 hr is much 
smaller than for CVs, as the lack of the enormous influx of newly-formed systems boosting the CV 
distribution there would suggest. Despite this lack the expected intrinsic period distribution (Fig. 4, 
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middle panel) derived from a typical evolutionary sequence (Fig. 4, upper panel) still predominantly 
populates the short period range, simply because the period evolves more slowly there (number 
density oc P^^)- However, these short-period systems are suppressed for any visibility function 
oc (— M2)" with a 1 (Fig. 4, lower panel). Such a visibility function corresponds for example 
to a flux-limited sample taken from a disc-like population. In addition, the detection probability 
function obtained in this way shows a pronounced peak at long orbital periods (here P 2± 8 hr), a 
feature consistent with the observed LMXB overpopulation at long orbital periods compared to CVs 
(see Fig. 3b). This peak is a consequence of the high mass transfer rate immediately after contact 
is reached in NLMXBs, since the systems are close to instability (M2 ~ Mi). The complexity of 
quantifying the relevant selection effects makes it difficult to decide if a ;^ 1 properly describes the 
NLMXB population. However, only if indeed a < 1 is an additional mechanism needed to account 
for the lack of short-period systems. One such mechanism is the evaporation of the secondary star 
by pulsar irradiation from a rapidly rotating neutron star spun-up by accretion (van den Heuvel &i 
van Paradijs 1988). Remarkably, this implicitly assumes that all LMXBs cross the gap, i.e. assumes 
the result we have demonstrated above. 



4. CONCLUSIONS 

We have shown that short-period neutron star low-mass X-ray binaries forming from helium- 
star supernovae without kick velocities must have secondaries with masses in the narrow range 
I.3M0 ^ M2 ^ 1.5M0. For standard strenghth magnetic braking the secondary star is already 
significantly nuclear-evolved when mass transfer begins, explaining why the resulting mass transfer 
rates are in many cases low enough for a substantial fraction of these systems to appear as soft 
X-ray transients even at short {P ^ 1 — 2 d) orbital periods (cf Paper 1). On the other hand, if the 
neutron star receives a strong kick velocity at birth, many NLMXBs would form with unevolved 
low-mass donors. The observed large fraction of SXTs among NLMXBs then forces us to conclude 
that kick velocities must be small compared to the pre-SN orbital velocity (^ 10%, i.e. ^ 50 km/s) 
for a large fraction of progenitor systems. This is consistent with a recent reevaluation of observed 
pulsar proper motions which suggest that the distribution of neutron star velocities at birth has a 
maximum at zero. Similarly, short-period NLMXBs forming from both initially thermally unstable 
systems (Kalogera & Webbink 1996a, b) and via the direct SN channel (Kalogera 1996b) would 
have a large fraction of (low-mass) unevolved secondaries, suggesting that neither of these channels 
contribute significantly to the short-period NLMXB population. 

Ignoring the uncertain survivors of thermally unstable mass transfer, the very special formation 
conditions for the case with negligible kick velocity also show that only a very restricted progenitor 
population (essentially ISM© + I.4M0 binaries with periods Pprog ~ 4 yr) can form NLMXBs, 
explaining their rarity in the Galaxy. We estimate a total formation rate ~ 2 X IQ-^ yr"^ The 
formation conditions are much more restricted than for cataclysmic variables. In particular short- 
period NLMXBs must all begin mass transfer at periods ^ 12 hr, in contrast to CVs, of which a 
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majority start mass transfer at periods ^ 1 — 2 hr. The resulting NLMXB period histogram has 
far fewer systems at short periods than the CV version, in agreement with observation. 

The smallness of the area of the Mne — M2 plane (Fig. 1) allowing short-period NLMXB formation 
is very striking. The fact that the resulting population has several properties in good agreement 
with observation is implicit confirmation that the assumed formation conditions are realistic. In 
particular it is clear that the neutron star mass at formation cannot be significantly larger than 
I.4M0: if it were, the allowed area would become much larger, sharply decreasing the predicted 
relative population of transient systems. However, we can say nothing about the lower limit on 
the formation mass, as systems with mi ^ 1.2 would not appear as NLMXBs (the allowed area 
in Fig. 1 would disappear). This result suggests another conclusion: with initial conditions 1.2 ^ 
mi ^ 1.4, 1.2 ^ 1712 ^ 1-5, evolution without mass loss would lead us to expect neutron star masses 
^ 2.4M0 in short-period systems, or those where direct estimates of m2 give a low mass, such as 
Cen X-4 (Shabaz et al. 1993). As there is no observational support for such masses, this suggests 
that a large fraction of the mass transferred in the NLMXB phase is lost from the binary. The 
most likely way for this to occur is through mass loss from the accretion discs in these systems (cf 
Begelman, McKee & Shields, 1983; Czerny & King, 1989). If the mass is lost at disc radii much 
greater than the size of the neutron star the central accretion rate can in principle be smaller than 
the mass transfer rate, so it may be simplistic to infer the latter from the observed X-ray flux. This 
in turn would mean that the presence of transient behaviour would pose a somewhat less stringent 
upper limit on the mass transfer rate than we inferred in Paper 1. 

This paper has discussed the formation of LMXBs containing neutron stars. The constraints on 
the formation of black-hole LMXBs appear to be weaker (cf, e.g., Romani 1994, 1996), as are the 
conditions for them to appear as transients (Paper 1). We shall investigate this problem in a future 
paper. 
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formation prior to publication and for very helpful discussions. We also thank the referee for useful 
comments. This work was supported by the UK Particle Physics and Astronomy Research Council 
through a Senior Fellowship (ARK) and a Rolling Grant for theoretical astrophysics to the Leicester 
Astronomy Group. 
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Fig. 1. — Constraints on NLMXB formation in the donor mass (M2) - He star mass (A/hc) plane in 
the absence of kick velocities. The straight line with positive slope represents the SN survival condition 
and is a lower limit for M2. The curve with negative slope represents another lower limit for M2 (or, 
equivalently, Mhg), the Galactic age constraint (for M2 > 0.3Mq; as magnetic braking does not operate 
in systems with fully convective secondaries we plot for AI2 < O.SMq the corresponding limit obtained by 
replacing ^mb by tcR.)- The horizontal line represents an upper (stability) limit for M2 (where 1.4Mo was 
assumed for the neutron star mass). These limits restrict the allowed range of NLMXB formation to the 
small triangular (shaded) area enclosed by these curves. Also shown is the lower limit for M2 (or A/ho) if the 
system's age is less than ^ms (dotted curve; / = 1 in eq [|3|), or less than 0.25iMS (dash-dotted; / = 0.25) 
at mass transfer turn-on. The allowed range for short-period NLMXBs with main-sequence donors is only 
the narrow-hatched area. 
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Fig. 2. — Same as Fig. 1, but for weak {£_ = 0.1, upper panel), moderate = 0.3, middle panel) and strong 
= 1.0, lower panel) kick velocities. For the ^ = 1 case the SN survival line was replaced by a line with 
P = const. = 0.4 (see Sect. 2.2). 
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Fig. 3. — Normalized cumulative period distribution for CVs and LMXBs in selected period ranges, data 
taken from Ritter & Kolb (1996). a: (upper panel) systems with orbital periods between 80 min and 2 d 
(249 CVs, 30 LMXBs); h: (lower panel) systems with orbital periods between ~ 3 h and 2 d (142 CVs, 27 
LMXBs). The arrows indicate the period where the difference between the distributions is largest. 
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Fig. 4. — Possible sliape of tlie NLMXB period distribution if all systems form with a secondary which is 

somewhat nuclear evolved and has a rather high mass. As a representatitve evolutionary sequence we use a 
calculation by Singer (Singer ct al. 1993, Ritter 1994) with an initial donor mass Ad2 = l-SM©, initial central 
hydrogen mass fraction 0.36 (i.e. with age = O.S^ms at turn-on), and a (constant) mass Mi = l.OM© for 
the compact object. The mass transfer rate M as a function of orbital period P for this sequence is shown 
in the upper panel. The sequence was originally used to represent CV evolution; a slightly different initial 
primary mass (I.4M0 instead of l.OM©) does not affect the main conclusions drawn in the text. The middle 
panel depicts the intrinsic discovery probability oc P/P vs log P, a quantity representing the intrinsic period 
distribution of a population of such binaries with similar initial configurations at mass transfer turn-on. In 
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